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^ Abstract 

P l This paper is concerned with a conservation law model of traffic flow on a network of 

roads, where each driver chooses his own departure time in order to minimize the sum of 
a departure cost and an arrival cost. The model includes various groups of drivers, with 
^3 different origins and destinations and having different cost functions. Under a natural 

set of assumptions, two main results are proved: (i) the existence of a globally optimal 
£^ solution, minimizing the sum of the costs to all drivers, and (ii) the existence of a Nash 

equilibrium solution, where no driver can lower his own cost by changing his departure 
I time or the route taken to reach destination. In the case of Nash solutions, all departure 

rates are uniformly bounded and have compact support. 
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1 Introduction 



Consider a model of traffic flow where drivers travel on a network of roads. We denote by 
Ai, . . . , A m the nodes of the network, and by 7^' the arc connecting Ai with Aj. Following the 
classical papers |13|. I14j. along each arc the flow of traffic will be modeled by the conservation 
law 

Pt + \pvij{p)] x = 0. (1.1) 

Here t is time and x £ [0, Lij] is the space variable along the arc 7^. By p = p(t, x) we denote 
the traffic density, while the map p h-> Vij(p) is the speed of cars as function of the density, 
along the arc 7y. We assume that v%j is a continuous, nonincreasing function. If Uy(0) > we 
say that the arc 7^ is viable. It is quite possible that two nodes i,j are not directly linked by 
a road. This situation can be easily modeled by taking = 0, so that the arc is not viable. 



The conservation laws (1.1) are supplemented by suitable boundary conditions at points of 



junctions, which will be discussed later. 

We consider n groups of drivers traveling on the network. Different groups are distinguished 
by the locations of departure and arrival, and by their cost functions. For k £ {1, ... ,n}, let 
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Gk be the total number of drivers in the fc-th group. All these drivers depart from a node ^(fc) 
and arrive at a node A a ( k \, but can choose different paths to reach destination. Of course, we 
assume that there exists at least one chain of viable arcs 

r = ( 7 i(0),i(l) ' 7 i(l),i(2) ' • ' ■ ' ( L2 ) 

with i(0) = d{k) and i{y) = a(k), connecting the departure node A d ^ with the arrival node 
A a nA. We shall denote by 

v = {r l7 r 2 , ... ,v N ] 

the set of all viable chains (i.e. concatenations of viable arcs) which do not contain any closed 
loop. Since there are m nodes, and each chain can visit each of them at most once, the 
cardinality of V is bounded by (m + 1)!. For a given k E {1, . . . , n}, let V k C V be the set of 
all viable paths for the fc-drivers, connecting Amq with A a ^ k y 

By U ktP (t) we denote the total number of drivers of the k-th group, traveling along the viable 
path r p , who have started their journey before time t. 

Definition 1. A departure distribution function t i— > U ktP (t) is a bounded, nondecreasing, 
left- continuous function, such that 

U kp (-oo) = lim U kp (t) = 0. 

Given group sizes G\, . . . , G n > 0, we say that a set of departure distribution functions {Uk, p } 
is admissible if it satisfies the constraints 

^£4, P (+oo) = G k k = l,...,n. (1.3) 

v 



Since G k is the total number of drivers in the fc-th group, the admissibility condition (1.3) 
means that, sooner or later, every driver of each group has to depart. If the function U k , p is 
absolutely continuous, its derivative will be denoted by 

Clearly, u kjP measures the rate of departures of fc-drivers traveling along T p . 

The overall traffic pattern can be determined by (i) the departure distribution functions U k ^ p (-), 



(ii) the conservation laws (1.1), and (iii) a suitable set of conditions at junctions, specifying 



the priorities assigned to drivers that wish to enter the same road. 

In this paper we consider the simplest type of condition at junctions, where a separate queue 
can form at the entrance of each road. Drivers arriving at the node Ai from all incoming 
roads, and who want to travel along the arc jij, join a queue at the entrance of this outgoing 
arc. Their place in the queue is determined by the time at which they arrive at Ai, first-come 
first-serve. Some additional care is needed to handle the case where different groups of drivers 
depart from the same node. Indeed, if a positive amount of drivers initiate their journey 
exactly at the same time, some additional information is needed to determine their place in 
the queue. This can be achieved in terms of the prioritizing functions introduced in [3]. 
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As in [21 El H] , we consider a set of departure costs </?&(•), and arrival costs ipk(') for the various 
drivers. Namely, a driver of the fc-th group departing at time T d and arriving at destination 
at time r a will incur in the total cost 



In this framework, the concepts of globally optimal solution and of Nash equilibrium solution 
considered in [21 [3] can be extended to traffic flows on a network of roads. 

Definition 2. An admissible family of departure distributions is globally optimal if 

it minimizes the sum of the total costs of all drivers. 

Definition 3. An admissible family {Uk tP } of departure distributions is a Nash equilibrium 

solution if no driver of any group can lower his own total cost by changing departure time or 
switching to a different path to reach destination. 

The main goal of this paper is to prove the existence of a globally optimal solution and of a 
Nash equilibrium solution, under natural assumptions on the costs and on the flux functions. 

In the case of a single group of drivers traveling on a single road, the existence and uniqueness 
of such solutions were proved in [2]. We highlight the main features of the present analysis. 

• Following the direct method of the Calculus of Variations, a globally optimal solution 
is constructed by taking the limit of a minimizing sequence {Uj?p} v >i of admissible 
departure distributions. The existence of the limit is guaranteed by the "tightness" of 
the sequence of approximating measures. Namely, for each e > there exists T > 
(independent of n) such that the total amount of drivers departing at times t ^ [— T, T] 
is less than e. 

• Toward the existence of a Nash equilibrium, the results in [21 [3] used the assumption tfj', > 
0, meaning that the arrival cost functions are nondecreasing. We now strengthen this 
assumption to ip' k > 0, so that the arrival costs are strictly increasing. This apparently 
minor change in the hypotheses has an important consequence. Namely, it allows us to 
prove a crucial a priori bound on all departure rates, in any Nash equilibrium solution. 

• The proofs in [21 [3] relied on a monotonicity argument. Indeed, the departure distribu- 
tion U(-) for a Nash equilibrium was obtained as the (unique) pointwise supremum of 
a family of admissible distributions, satisfying an additional constraint. On the other 
hand, the present existence result is proved by a fixed point argument. By its nature, 
this topological technique cannot yield information about uniqueness or continuous de- 
pendence of the Nash equilibrium. 

The paper is organized as follows. In Section 2 we describe more carefully the traffic flow 
model, explaining how to compute the admissible solutions using the Lax-Hopf formula. In 
Section 3 we prove the existence of a globally optimal solution, while Section 4 is devoted to 
the existence of a Nash equilibrium solution. 

For the modeling of traffic flow we refer to [H El H~3l E] . Traffic flow on networks has been the 
topic of an extensive literature, see for example El El [10] and references therein. A different 
type of optimization problems for traffic flow was considered in |llj . 
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2 Analysis of the traffic flow model 



In our model, x G [0, Lij] is the space variable, describing a point along the arc 7^. Here 
measures the length of this arc. The basic assumptions on the flux functions Fij(p) = pvij(p) 
and on the cost functions (pk,ipk are as follows. 

(Al) For every viable arc 7y, the flux function p \-t Fij(p) = pvij(p) is continuous, concave 
down, and non-negative on some interval [O,/?^], with Fij(0) = Fijiji^) = 0. We shall 
denote by p*j G]0,^ [ the unique value such that 

F(P*ij) = F™ X = max ^(p), F^p) > for a.e. p€[0,p£]. (2.1) 

pe[a,pij] 

(A2) For every k S {l,...,n} the cost functions tpk,ipk are continuously differentiable and 
satisfy 

lim (<p k (t)+Mt)) = +00. (2.2) 




Figure 1: Left: the function p 1— » Fij(p) = pvij(p) describing the flux of cars. Middle: the function 
gij, implicitly defined by gij(pvij(p)) = p. Right: the Legendre transform g*j. 



Remark 2.1. By (Al), the flux function u = Fij(p) is continuous, concave and strictly 
increasing on the interval [0,p*-]. Therefore, it has a continuous inverse: p = gij(u). As shown 
in fig. [TJ the function u h-> gij{u) is convex and maps the interval [0, F™ ax ] onto [0,^*^]. 

Remark 2.2. According to (A2), the cost for early departure is strictly decreasing in time, 
while cost for late arrival is strictly increasing. The assumption that these costs tend to infinity 
as t — > ±00 coincides with common sense and guarantees that in an equilibrium solution the 
departure rates are compactly supported. 

Remark 2.3. In the engineering literature (see for example [9]) it is common to define the 
travel cost as the sum of the travel time plus a penalty if the arrival time does not coincide 
with the target time Ta ■ 

D(t) + H>(t + D(t)-T A ). (2.3) 

Here D{t) denotes the total duration of the trip for a driver departing at time t, while ^ is a 
penalty function. Calling r a {t) = t + D(t) the arrival time of driver departing at t, the cost 
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function (2.3) can be recast in the form (|L5j). Indeed, 

D(t) + V(t + D(t)-T A ) = -t + t + D(t) + 9[t + D(t)-T A ) 

= -t + T a {t) + y(T a (t)-T A ) 
= ^)+V(T a (t)), 

where (p(t) = —t, iP(t) = r + \&(t — T A ), In order that the assumptions (A2) be satisfied, it 
suffices to require that the function \P be continuously differentiate and 



* > 0, 



1, 



lim - t) 

\t\— >oo 



+ OO . 



(2.4) 



2.1 Traffic flow with an absolutely continuous departure distribution. 

We now describe more in detail how the traffic flow on the entire network can be uniquely 
determined, given the departure distributions p . As a first step, we consider the absolutely 



continuous case, so that (1.4) holds. 

Along any arc 7y, the traffic density pij satisfies a boundary value problem of the form 
(d t pij(t, x) t + d x F ij (p ij (t, x)) =0 (t, x) G R x [0, Ly] , 

jFiiOM*, °)) = 1 eM > 

where u~At) describes the incoming flux at x = 0. Following the approach in [21 13], we switch 



the roles of the variables t, x, replacing the above boundary value problem (2.5) with a Cauchy 
problem for the conservation law describing the flux ity = Fij(pij)' 

(d x Uij{t, x)+d t gij(uij(t, x)) =0 (t, x) G K x [0, Ly] , 

\ Uij (t, 0) = u^it) tel. 



Here g%j : [0, F^ ax ] h-» [O,/^-] is the inverse of the function Fij, as in Remark 2.1 Consider 
the integrated functions 

Uij(t, x) = \ Uij(s,x)ds, UfAt) = / vZ(s)ds. 



00 



Then Uij(t, x) provides a solution to the Hamilton- Jacobi equation 

(d x Uij(t, x)+g ij (dtU ij (t, x)) =0, 
{UyfaO) = Ur.(t). 

The viscosity solution to the above Cauchy problem is given by the Lax-Hopf formula [7] 

Uij(t, x) = mm I U7j(r) +xg* ij (^p) \ , x G [0, Ly] . (2.8) 
Here <?|j- denotes the Legendre transform of gij, namely 

9ij (p) = r max , [Pu-gij(u)). (2.9) 
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Figure 2: A generic node of the network. 



For any given node Ai, in general there will be several incoming arcs 7^, I S I(i) and several 
outgoing arcs 7^, j € 0(i) (see fig. [2|. To determine the flux at the entrance of the arc 7^, 
one needs to know how many drivers, after reaching the node A4, actually want to take the 
road jij. For this purpose, we need to introduce the distribution functions 



• U kpi j(t) = total number of drivers of the k-th group, traveling along the path T p , that 
have entered the arc 7™ (possibly joining a queue at the entrance) within time t. 

• U^ pi -{t) = total number of drivers of the k-th group, traveling along the path T p , that 
have exited from the arc 7^ (reaching the node Aj) within time t. 



Of course, U kp - = if the path T p does not contain the arc 7^. 

The entire flux along the network is entirely determined by the functions Uj^ - . To recursively 
compute these functions, observe that the total number of drivers who have entered the arc 
jij (possibly joining a queue) within time t is determined by 

W= EE E E E (2-io) 

eex(i) k yn,7ijEr p d{k)=i^j& p 

Notice that the first sum accounts for the drivers that transit through the node Ai, while the 
second sum accounts for the drivers that initiate their journey from the node A}. Using the 
solution formula (2.8) with x = Lij we obtain 

U$(t) =U l3 {t,L l0 ) = m m{l^(T) + L^(^)} . (2.11) 

Notice that the function Uf- is nondecreasing and Lipschitz continuous. Indeed, the rate at 
which drivers arrive at the end of the arc jij is computed by 

= jUtjit) G [0,f$H fora.e. t. 

Among the drivers who reach the end of the arc 7y within time t, we still need to compute 
how many belong to the various groups. 

In the case where all departure rates Uk tP with d(k) = i are absolutely continuous, the functions 
J7i ■■ can be computed as follows. For a.e. t there exists a unique time T enter (t) such that 



U- Renter {t)) = jj+ {t y 
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A driver entering the arc 7y at time T enter {t) thus reaches the end of the arc at time t. The 
first-come first-serve assumption now implies 



= gf^jf u ' eJ(/ ' ) - --"--^r,- (-'-12) 



In other words, among all drivers who exit from the arc 7^ at time t, the percentage of (k,p)- 
drivers must be equal to the percentage of (A;,p)-drivers that enter 7« at time T enter (t). In 
turn, the arrival distributions are computed by 

Ui, l3 (t) = f u+#(f)cft. (2.14) 

J —00 

Remark 2.4. The above equations can be solved iteratively in time. Namely, let 

. U 



'mm 



-'11 

mm 



%(0) 

be the minimum time needed to travel along any viable arc of the network. Given the departure 
rates Uk,p, if the functions Ujz, .-(i) are known for all i < r, by the above equations (2.10)-(2.14) 
one can uniquely determine the values of U^ p -(t) also for t < r + A m « n . 

Let Gk, P be the total number of drivers of the k-ih group who travel along the path T p . 
The admissibility condition implies Gf-i + • • • + Gk,N = G^. We use the Lagrangian variable 
/3 6 [0, Gk p ] to label a particular driver in the subgroup Gk, P of /c-drivers traveling along the 
path F p . The departure and arrival time of this driver will be denoted by T^ p {j3) and t% (/3), 

respectively. Let U kp art (t) = Uk. p {t) denote the amount of drivers of the subgroup Qk,p who 
have departed before time t. Similarly, let U k r p we (t) be the amount of (&;,p)-drivers who have 
arrived at destination before time t. For a.e. (3 £ [0, Gk, p ] we then have 

7&G8) = inf{r; U**"*® > fi) , < p (/3) = inf{r; U%™ e (t) > ft}. (2.15) 

With this notation, the definition of globally optimal and of Nash equilibrium solution can be 
more precisely formulated. 



Definition 2'. An admissible family of departure distributions {Uk tP } is a globally optimal 
solution if it provides a global minimum to the functional 

J = E f Gk '\^k{ri p m+Mrt p m)dP- (2.16) 
k,p J ° 

Definition 3'. An admissible family of departure distributions {Uk tP } is a Nash equilibrium 
solution if there exist constants c\, . . . , c n such that: 

(i) For almost every (3 £ [0,Gk, p ] one has 

<Pk(Tl p m+Mr£ >p m = (2.17) 
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(ii) For all t £ M, there holds 



(2.18) 



Here Aj. jP (r) is the arrival time of a driver that starts at time r from the node and 
reaches the node A a ^ traveling along the path T p . 

In other words, condition (i) states that all £;-drivers bear the same cost c&. Condition (ii) 
means that, regardless of the starting time x, no fc-driver can achieve a cost < c^. 



2.2 Traffic flow with general departure distribution 

If all departure distributions Uk tP are absolutely continuous, the previous analysis shows that 
the first-come first-serve assumption on the queues completely determines the traffic pattern. 
This is no longer true if a positive amount of drivers of different groups initiate their journey 
exactly at the same time. For example, assume that i = d(k) = d(k') is the departure node for 
both fc-drivers and fc'-drivers. Let jij be the first arc in the paths T p and V p /, and assume that, 
at the instant to, a positive amount of drivers in the subgroups Gk,p and Gk',p' initiate their 
journey. Since all these drivers join the queue at the entrance of the arc 7^ at the same time, 
additional information must be provided to determine their relative position in the queue. For 
this purpose, we follow the approach introduced in [3] . 

Given an arc 7^, call 

Gij = |J Gk, P (2.19) 
d(k)=i, 7^er p 

the family of all drivers that initiate their journey from the node A4, traveling along 7^ as 
first leg of their journey. The total number of these drivers is 

Gij = ^ ^ Gk, p ■ 
d(fc)=i,7ijer p 

The cumulative departure distribution u^ epart : M. 1— > [0, Gij] is defined as 

uf; part (t) = J2 u ^)- ( 2 - 2 °) 

In other words, uf^ paTt {t) is the total number of drivers in the family Gij that depart before 
time t. The relative position of these drivers in the queue at the entrance of the arc 7^ will 
be determined by a set of prioritizing functions. 

Definition 4. A set of prioritizing functions for departures on the arc jij is a family of 
nondecreasing maps 

Bk, P '■ [0, Gij] H> [0, Gk, p ] , Gk, P Q Gij , (2-21) 

satisfying 

B k , p ((3) = (3 for all 6 [0,G lJ ], (2.22) 



s 



B kiP (U% part (t)) = U k;P (t) fora.e. t€R. (2.23) 



Otherwise stated, among the first /3 drivers that depart and choose 7^ as the first arc of 
their journey, £>fc jP (/3) counts how many belong to the subgroup Gk,p- If all functions Uk, P 



are continuous, then the conditions ( 2.22 )-( 2.23) uniquely determine the prioritizing functions 
B k>p . In this case, there is actually no need to introduce this concept. On the other hand, 
if two or more functions U k , p have an upward jump at a time to, then different prioritizing 
functions are possible. 

We now show that, by choosing one set of prioritizing functions, the traffic flow on the entire 
network is entirely determined. Indeed, for all k,p,i,j the values U^ pi -{t), determining how 
many drivers of the subgroup Gk,p reach the end of the arc 7^ within time t, are computed as 
follows. 

CASE 1: The only drivers traveling on the arc 7^ are those who started their journey from 



the node A{. In this case, as in (2.11) the total number of drivers arriving at the end of the 
arc jtj within time t is 



E = m . in {^w + ^^(V)} ' (2,24) 



k,p 

Given the prioritizing functions Bk, p , the values U^ pi -{t) are immediately obtained by the 
formula 

= • ( 2 - 25 ) 

CASE 2: The arc 7y is also traveled by drivers who transit through the node Ai, departing 
from other nodes. In this case, drivers originating from Ai have to merge with drivers in 
transit, coming from other nodes. The cumulative distribution function, accounting for the 
total number of drivers that have entered the arc 7^ within time t is 

U~ 3 {t) = U^ a ' rt (t) + Ul™™\t) = £ U^(t)+ £ U^{t). (2.26) 

As before, the total number of drivers who have exited from the arc 7^ before time t is given 
by (2.11). To determine how many of these drivers belong to each subgroup G k , P , we proceed 



as follows. 

Consider the driver who exits from the arc 7^ at time t. This driver will have entered the arc 
jij at an earlier time r = T enter (t), such that 

Vim ._Ur.( 3 ) < U+(t) < Vn^Ur^s). (2.27) 

Notice that T enter (t) is uniquely determined, for all but countably many times t. For all 
subgroups of drivers in transit, the first-come first-serve priority assumption implies 

U^jit) = Uj^^m (2.28) 



Here the right hand side of (2.28) is uniquely determined, because the distribution function 
^kpij m (2-26) is Lipschitz continuous. 
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On the other hand, for drivers who initiate their journey at A{, the distribution function U kp - 
can have a jump at r enter (i), in which case the formula (2.28) is not meaningful. Given a set 



of prioritizing functions B k>p , the distribution functions can be determined by setting 

P = U±(t)-U% ansit (T enter (t)), (2.29) 

K,H = B >*M ■ (2-30) 

Indeed, U^(t) is the total number of drivers that reach the end of the arc 7^ before time 
t, while p counts how many of these drivers start their journey from the node A^. In turn, 
B k<p {J5) determines how many belong to the subgroup Q k;P . 



3 Globally optimal solutions 

In this section we establish the existence of a globally optimal solution. The proof follows the 
direct method of the Calculus of Variations, constructing a minimizing sequence of solutions 
and showing that a subsequence converges to the optimal one. 

Theorem 1 (existence of a globally optimal solution). Let the flux functions and the 
cost functions ip k ,if) k satisfy the assumptions (A1)-(A2). Then, for any n-tuple (G\, . . . ,G n ) 
of nonnegative numbers, there exists an admissible set of departure distributions U kp and 
prioritizing functions B kp which yield a globally optimal solution of the traffic flow problem. 



Proof. 1. By (A2), all functions ip k + ip k are bounded below. By possibly adding a constant, 
it is not restrictive to assume that ip k (t)+ijj k (t) > for every time t. Calling mo the infimum of 
all total costs in (2.16), taken among all admissible departure distributions {U k ^ p }, this implies 
mo > 0. In addition, it is clearly not restrictive to assume Gk > for all k £ {1, . . . , n}. 

Recalling Definitions 1 and 4, consider a minimizing sequence of departure distributions U% 
and prioritizing functions B kp : [0, G^] t-t [0,G kp \. Here 

G t, P = U k, P (+°°)- 
By choosing a subsequence, we can assume 

lim G u k = G kp with y^G kp = G k . (3.1) 

V 

Moreover, by Helly's compactness theorem we can assume that, as v — > 00, one has the 
pointwise convergence 

U^ p {t) -»• U k>p (t) for a.e. t el, (3.2) 

while Ascoli's theorem yields the uniform convergence 

Bl p ((3) -> B k , p {P) peR. (3.3) 



In (3.3), we have extended the prioritizing functions to the entire real line by setting 

[0 if /3<0, 
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In the remainder of the proof we show that the set of departure distributions Uk, P together 
with the prioritizing functions B k>p yield a globally optimal solution. 



2. Let e > be given. We claim that there exists a large enough constant T, independent of 
v such that 

' £c/£„(i) < e for t<-T, 

p 



> G k -e 



(3.4) 



for t > T , 



for all v sufficiently large. Indeed, by (A2) there exists T such that 

m + 1 



for \t\>T, k 6 {!,..., n} 



If any one of the two conditions in (3.4) fails, then the total cost would be > mo + 1. Since 
by assumption as v — > oo the total cost approaches the infimum mo, this proves our claim. 



By (3.4) it follows that the limit functions satisfy 

U kjP (-oo) = 0, U kjP (+oo) = G k , p 

In particular, the limit departure distribution is admissible. 



3. For G [0, G v k p ] let 

P ^ and ^ r^((3) (3.5) 

describe the departure and arrival time of the /3-driver, in the subgroup Qk, P - We claim that, 
by possibly extracting a further subsequence, one has the pointwise convergence 

t&(P) < P (P)^ < P (P) fora.ejG[0,GJ. (3.6) 



Indeed, the maps in (3.5) are nondecreasing. Moreover, for any e > 0, by (3.4) these maps 
are uniformly bounded when restricted to the subinterval [e, Gk, P — s], uniformly w.r.t. v. 
By Helly's compactness theorem, we can find a subsequence that converges pointwise on 
[e, Gk, P — s]. Since e > was arbitrary, a standard argument proves our claim. 

4. It remains to prove that the limit departure distribution is optimal. Recall that, without 
loss of generality, we are assuming ipk(t) + ifrkit) — f° r an k,t. The total cost J determined 
by the set of departure distributions {Uk p } and prioritizing functions {Bk p } can be estimated 
by 



J = E/ ^ (Mri P (P)) + M^, P (P)))d(3 
k,p Jo 

(■>-) 

= ™P E/ ^ £ (^« P (/ ? ))+^« p (/3)))d/3 
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Fix e > 0. By (3.4) all departures and arrivals of /3-drivers with /3 G [e, Gk lP — e] take place in 
a uniformly bounded interval of time, say [—T,T']. On this interval, all functions (pk,ipk are 
uniformly continuous. Hence the pointwise convergence (3.6) yields 



£/ '"" (vk{Tt P {i3))+i>M, P m)dP 



h::p 





= m . 



Together with (3.7), this implies J < mo, completing the proof. 



□ 



Remark 3.1. The above theorem remains valid if in assumption (A2) we only require ip' k < 0, 



Remark 3.2. A natural conjecture is that, in a globally optimal solution, all departure rates 
Uk tP = ^Uk tP are uniformly bounded and have compact support. Hence the corresponding 
solution should be uniquely determined without need of prioritizing functions. In the case 
of one group of drivers traveling on a single road, this fact was proved in [2J. In the general 
case, a proof of the above conjecture will likely require a more detailed study of the globally 
optimal solution, establishing necessary conditions for optimality. 



4 Nash equilibria 

In this section we prove the existence of a Nash equilibrium solution for traffic flow on a 
network. For our model, it turns out that in a Nash equilibrium all departure rates must be 
uniformly bounded and have compact support. As a consequence, the corresponding solution 
is uniquely determined without need of prioritizing functions. 

Theorem 2 (existence of a Nash equilibrium) . Let the flux functions and the cost 
functions (pkirfk satisfy the assumptions (Al)-(A2). 

(i) For any n-tuple {G\, . . . , G n ) of nonnegative numbers, there exists at least one admissible 
family of departure rates {u* k p } which yields a Nash equilibrium solution. 

(ii) In every Nash equilibrium solution, all departure rates are uniformly bounded and have 
compact support. 

Before proving the theorem, we establish a "modulus of continuity" for the exit time. Namely, 
for drivers traveling along a given path r„, the arrival time r p (t) is a uniformly continuous 
function of the departure time t. In the following, we first consider a single arc 7^ with flux 
function Fij(-) satisfying (Al). As in figure pi gij denotes the inverse function while is the 
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Legendre transform. For a driver who enters the arc 7^ at time t (possibly joining a queue), 
we denote by Tij(t) his exit time. 



Lemma 4.1. Given constants G, M > 0, there exists a continuous function <f>. 



depending on the flux function Fij in (2.1) and on M,G, such that (f>ij(0) = and moreover 
the following holds. Let U~A-) be a Lipschitz continuous departure distribution, such that 
< Uij(t) = ^U~-{t) < M for a.e. t and J Uij(t) dt < G. Then the exit time %■(•) satisfies 

Tij(h) — nj(ti) < <pij(h — ti) whenever ti < t 2 ■ (4.1) 



Proof. 1. Let Ly be the length of the arc 7y. Following [2], consider the function 

hijia) = ' ( 4 - 2 ) 

Since the Legendre transform g*- is convex, the function hij is concave. Moreover, calling 
fj,ij = Lij/vij(0) the minimum time needed to drive across the arc 7^ (= length of the arc 
divided by the maximum speed), one has 



hij(s) =0 if s > -mj , 
hij(s) < if s < —fiij ■ 



(4.3) 



Using the solution formula (2.24), the amount of drivers that exit from the arc 7^' before time 
t is computed by (see Fig. [3| 



mm 

s 



{u^-h^s-r)). 



(4.4) 



In other words, U^(t) is the amount by which one can shift upward the graph of h(- — r), 
before hitting the graph of U^(-). In turn, the exit time of a driver who enters the arc 7y at 















IJ 

/ y 


t T.(t) 

U 



Figure 3: A geometric construction of the exit distribution t M> U^(t) from the entrance distribution 
U^j , using the formula (4.4). 



time t is given by 

nj(t) = (t + Hij) V inf {r; V~ (a) > Vr. (t) + (s - r) for all s < i} , (4.5) 
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where we used the notation a V b = max{a, b}. 

2. Let t\ < t 2 be given. Two cases can be considered. 

CASE 1 (Fig. left). The minimum 

Utjiratb)) = min {ur.{ a )-hij{8-ni{h))}. (4.6) 
is attained at a point s G [ii, ^2]- In this case we have 

7"y(*l) > h + Mij i 
Tij(t 2 ) <t2 + inf|r; Ms>hij(s-T) for all s € [-M(i 2 - *i) , 0]} . 

Hence 

Tij(t 2 ) -Tijih) < ^-(ia-ii) (4.7) 

with 

4j(0 = e + mf {r ; Ms > fcy(a - r) for all s G [— M£ , 0]} - . (4.8) 



CASE 2 (Fig. [4} right). The minimum in (|4.6l) is attained at a point s < t\. In this case we 
have 

Ur.itJ + hijis-njih)) < Ur.(s) = u^h) + h^s - Tij (t 2 )) . 
Since hij is concave and Tij{t\) —t\> fiij, Tij(t 2 ) — t 2 > fJ-ij, one has 

> hij (Tij(ti) - Tij(t 2 ) - (iij) > hij(s - Tij(t 2 )) - hij(s - Tij(h)) 
In this case, we conclude 

Tij(t 2 ) - njitx) < fyh-h), (4.9) 
where the continuous function (jr- is implicitly defined by 



3. The two functions 4>\p (fr^ defined at (4.8), (4.10) are both continuous and vanish at the 



origin. Defining 

<M0 = max 4j(0}, (4.11) 

the conclusion of the lemma is satisfied. □ 
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t, 8 t 2 Tj/tt) 




Figure 4: The two cases considered in the proof of Lemma 4.1 



The next lemma extends the result of Lemma 4.1 to a general path T p . Here r p (t) denotes 
the arrival time of a driver starting at time t and traveling along T p . 



Lemma 4.2. Let all departure rates Uk, p (t) be uniformly bounded, so that 

Uk, P (t) < Mq for all k,p,t. 



(4.12) 



Then, for any viable path T p , there exists continuous function <p p : M + i— > K + such that 

4> p (0) = 0, ( T p(h) — Tp(h)) < <j> P (h — ti) whenever ti<t 2 - (4.13) 



Proof. The assumption (4.12), together with the fact that the flux through each arc jij 
cannot be greater than the maximum flux F™ ax , implies that the maximum incoming flux 
through each arc is bounded by some constant M. For each arc 7^, let be the modulus of 



continuous dependence constructed in Lemma 4.1 If T is the path in (1.2), obtained as the 



concatenation of the arcs JM-i)M£), £ = 1, • ■ ■ , v, it now suffices to define the function <f> p as 
the composition of the corresponding functions (f>M—i),i(£)i nam ely 



5 i(l),i(2) <W(0),i(l) 



□ 



Our final lemma shows that the arrival times, and hence the cost functions, depend continu- 
ously on the departure rates. 

Lemma 4.3. Consider a sequence of departure rates u v = (11% ) which are uniformly bounded 
and supported inside a common interval I = [—T,T], namely 



< u% >p (t) < Mo for all t€l, uj^t) = if t <£ I , 

for every v > 1. Assume that for all k,p one has the weak convergence 

u k,p 



II 



k,p • 



(4.14) 
(4.15) 



For each viable path T q , call T^{t), r*{t) the corresponding arrival times of a driver who departs 
at time t and travels along T q . Then, as v — ^ 00, one has the uniform convergence 



T q ~ r q \\C([-T,T]) 



0. 



(4.16) 
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Proof. 1. We first consider the case of a single arc 7 = 7^. Assume that the departure rates 
u v satisfy 

«"(*) G [0,M ] if t G [— T, T] , 

u"(t) = if t£ [-7,7], 

and converge weakly: u 1 " — 1 u. Define the integrated functions 



U u {t) = f u u (s)ds, U*(t) = I u{s) 

J— 00 J — 00 



(is . 



Our assumptions imply the convergence U v (t) — > U(t), uniformly for t £ I. In turn, this 
implies that the exit distributions 

u»(t) = mm lirw+Li, > m<) = m T in {^)+^4(*xr)} > 

U (4.17) 



satisfy the uniform convergence E/+ — > 17+ . Indeed, by (4.17) we have 

\\U1 - £7+|| L cx> < || 17" - ?7|| L oc -+ as ^00. (4.18) 



2. Given e > choose v large enough so that 
\U v {t)-XJ{t)\ < e 



for all v > v, , t G 



(4.19) 



Fix a time t G [— 7, 7]. To estimate the difference \r u (t) — r(t)\ between the corresponding 
arrival times, consider the modified departure distribution 



17(a) if s<t, 

t/»(s) = { U(t) + (l + M )(s-t) if se[t,T], 
U(t) + (1 + M )(T-t) if s>T. 



(4.20) 



Call r^( t) th e arrival time of a driver departing at time i, relative to the distribution {/". By 



Lemma 



4.1 



the function t \-¥ T*(t) satisfies a uniform modulus of continuity say 
r tt (t 2 ) -r tt (ti) < 0(t 2 -ti) foralUi<t 2 , 
for some continuous function with (f)(0) = 0. In particular, 

r\t + e) -T*(t) < 4>{e). 

Observing that 

U\t + e) - U\s) > U u (t + e)-U u {s) for all v > v , s<t + e 



(4.21) 



and using the representation formula (4.5) with hij given by (4.2), we obtain 
r*(t + e) 



= (t + e + Wi ) V infjr; U* (t + e) - (s) < 

> (t + e + Wi ) V infjr; U v {t + e) - U v {s) < L tJ g\ 
= T»{t + e). 



-S+T 



for all s < t + ej 
for all s < i + e| 
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Therefore 

T v (t) < T u (t + e) < r'(i + e) < r\t) + <f>(e) = r{t) + cp{e) forallzy>z>. (4.22) 

Switching the roles of U, U u and hence defining 

( U u (s) if s<t, 

V*(a) = I U v (t) + (l + M )(s-t) ifse[t,r], 
[ U u {t) + (l + M )(T-t) if s>T, 

we obtain the converse inequality 

r{t) < r v {t) + <f>(e) for all v > v . (4.23) 



3. For each arc 7^, the inequalities (4.22)-(4.23) show that the functions t 1— )• converge 
uniformly to the corresponding functions 1 1— > Tij(t). Given any path r g , the uniform conver- 
gence Tg — > t* is now obtained by a straightforward argument, using induction on the number 
of arcs contained in T q . □ 

We are now ready to prove the main result of this paper. 
Proof of Theorem 2. 

1. We claim that there exists a time interval / = [—To, To] so large that, in any Nash equilib- 
rium, no driver will depart or arrive at a time t ^ I. Indeed, given the n-tuple (G\, . . . , G n ), 
the travel time along any viable path T p = (7i(o),i(i)> ■ • • , 7i(u-i),i(u)) is a priori bounded by 

VT = t { + „ L V^ m , } ■ (4.24) 

Here and in the sequel, we call 

G = G 1 + ... + G n (4.25) 

the total number of drivers. Notice that, in each summand on the right hand side of ( |4.24[ ), 
the first term is an upper bound for the time spent waiting in the queue (total number of 
drivers divided by the maximum flux) while the second term is an upper bound on the actual 
travel time (length divided by the minimum speed). Let 

T max = max T max 
p p 

be an upper bound on the travel time along all viable paths. In view of assumption (A2), 
there exists To large enough such that 

mm{Mt) + Mt)} > max{M0)+MT max )} for alU £ I = [-T , T ] . (4.26) 

k k 

Therefore, in a Nash equilibrium no driver will depart or arrive outside [—To, To]. Otherwise, 
he would achieve a strictly lower cost by departing at time t = 0. 

2. Let F max = maxF^ ra be an upper bound for the flux over all arcs. Call 

Vmax = max max \cp' k (t)\ , ij/ min = mm mm i// k (t). 

Kk<n t£l l<k<n t£l 
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Observe that V'mm > 0> because of the assumption (A2). We claim that, in a Nash equilibrium, 
all departure rates Uk p must satisfy the priori bound 

/.\ ^ ■ fmax ' F ma x c j. I a T7\ 

u k,p{t) < k = 77 for a.e. t. (4.27) 

Train 

Indeed, consider drivers of the fc-th family traveling along the path T p . Let t\ < t 2 be any two 
departure times, and call r% < t 2 the corresponding arrival times. The total costs for these 
two drivers must be the same, hence 

ifih) + tp{n) = (p(t 2 )+ip{T 2 ). 

On the other hand, the upper bound on the flux implies 



'■t'Z 

max Jt\ 

Therefore 



1 f t2 

T2-n > — / u k , p {t)dt. 

Fmax JU 



Fornax Jt 



*2 



1 



(t 2 -h)<p' max > <p(h) - (pfa) = ^(t 2 )-^(ti) > fo-Ti)VC > %™ I u kiP (t)dt. 
We thus conclude 

Vmax-Fmax > _J_ ^ , f) ^ _ 



Cm *2-*l 



Since this bound is valid for every interval , ^2] Q I, the pointwise bound (4.27) must hold. 
Moreover, for t £ I we already know that v,k tP (t) = 0. The last statement of the Theorem is 
thus proved. 



3. Choose the time 

T = T + -, (4.28) 

K 

where G = Y, k Gk ■ 

Consider the family of admissible departure rates 



U = \ {u kyP )i< k < n) i< p < N ; u k , p : M h-> [0,4k] , Uk,p(t) = for t <£ [-T,T] , 
J Uk, p {t)dt = Gk for every k X. 



(4.29) 

It is understood that Uk iP = if the path T p does not connect Amja with A a n.y Notice that 
U is a closed convex subset of L-^M; W nxN ). 

For each fixed > 1, we consider a finite dimensional subset U u CU consisting of all u = (uk tP ) 
which are piecewise constant on time intervals of length T j v. Introducing the points 

t 

t\ = -T , -v<l<v , 
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we thus define 

U v = jit = (uk )P ) G ; 

every function Uk, p is constant on each subinterval I\ = ]tg_i, 

(4.30) 

Observe that every u £U U has the form 

u = {u k ,p), u k ,p(t) = u k ,p,e for all (4.31) 

for some constants Ufc^ G [0, 4k]. 

4. Given u = (iifc iP ) £ £7, let r 9 (i) be the arrival time of a driver starting at time t and traveling 
along the path T q . Clearly, this arrival time depends on the overall traffic conditions, hence 
on all functions Uk tP - If this driver belongs to the j-th family, his total cost is 

*&?(*) = w(0 + ifofo,(f)). 

We now observe that, for each v > 1, the domain U u is& finite dimensional, compact, convex 



subset of L 2 ([-T,T]; R nxN ). Moreover, by Lemma 4.3 the maps u H- are continuous 

from Z7j, into L 2 . Hence, by the theory of variational inequalities |12j . there exists a function 
u 1 " = (u v a „) G Z7 y which satisfies 

/T 
^(t) • " «£,(*)) tft > for all veU v . (4.32) 

i-'i ' T 



5. We now let v — >• oo. By the previous steps, there exists a sequence of piecewise constant 
functions u v = (u kp ) G £7^ such that (4.32) holds for every v > 1. Since all functions are 
uniformly bounded and supported inside the interval 7 = [— T, T] , by taking a subsequence we 
can assume the weak convergence 

K P ) ^ «„) (4-33) 

for some function u* = ) G Z7. We claim that the departure rates u* kp yield a Nash 
equilibrium solution. More precisely: 



(NE) Given any k,p, any t\ G Supp(u* kp ), ti G M and any path T g with the same initial and 

n, P w ^ ( 4 - 34 ) 



final nodes as T p , one has 



Indeed, (4.34) implies that no &:-driver can lower his own cost by switching to the time t2 or 
choosing the alternative path T q to reach destination. We recall that t is in the support of a 
function / G L 1 if and only if J t ^ \f(s) \ ds ^ for every e > 0. 



6. By Lemma 4.2, all maps t i— > t% (t) have the same modulus of continuity. Since these 
arrival times are uniformly bounded, we can apply Ascoli's compactness theorem. Choosing 
a subsequence and relabeling, as v — > oo we thus achieve the convergence 

T k,p(t) ~~ ^ T k,p(f) for all k,p, uniformly for t G [—T,T]. (4.35) 
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By the assumption (A2), the cost functions <pk{ m ), ipk( m 

converge to 



tions ^(-) = Vk(-) + Mr^(-)) 
t € [— T, T]. By Lemma 
t and following the path 



4.3 

pi 



are continuous. Therefore, the func- 
= ^fc(-) + ^jfc(^|, p (0)> uniformly for 
TfcJt) is indeed the arrival time of a /c-driver departing at time 
in the case where the departure rates of all drivers are given by 



4k 



k,p 



4k 



'k.p 



I v 

i* 



4k 



'k, P 



k, P 



1 1 

j-i J 



Figure 5: The three cases considered in the proof of Theorem 2. In Case 1 (top left) the average 
of the cost $fc iP on the interval 1\ is higher than on the interval Jj\ To obtain a contradiction with 
(4.32) we simply move some of the mass from I\ to JJ. In Case 2a (top right) one cannot increase the 
value of u v k p on the interval because of the constraint u < 4k. However, some mass can be moved 
from I" to the previous interval In Case 2b (bottom) there are several adjacent intervals where 

u\ = 4k. In this case, if 2J» is the first interval to the left of 7J where m£ < 4k, we argue that (i) 
tj, > — T, and (ii) the average of the cost &k,p on ij« is strictly less than on ij. In this last case, to 
obtain a contradiction with (4.32) we move some mass from 2J„ to ij. 



7. If (4.34) fails, then by continuity there exists <5 > such that 

$lj,(t) > + 2(5 whenever |i - ti| < 2(5, |t' - i 2 | < 2«5 . 

By uniform convergence, for all v large enough we have 



whenever \t — t%\ < 2(5, Ii — ^ 2 1 < 2(5 . 



(4.36) 



(4.37) 



Observe that it is not restrictive to assume that t% £ [— Tq,Tq]. Indeed, if (4.34) fails for some 



*2 £ RZb, T ], then Q4.26[ ) implies 

and we can simply replace ti with zero. 



The weak convergence (4.33), together with the assumption on the support of the function 
u* k , now implies 



lim 

v— >oa 



h+8 



l k,p 



[t)dt 



h+s 



u* k Jt)dt > 0. 
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Therefore, for every v sufficiently large we can find two intervals 

I? = c fa- 5, h + 6], 

with V. > t 2 and u v k p {t) > for t G I? . 



(4.38) 



8. We now derive a contradiction, showing that, for v sufficiently large, the departure rates 
u kp do not satisfy the variational inequality (4.32). Two possibilities can arise (see figJHb. 



CASE 1: u kq • < 4k. In this case we define a new set of departure rates v 6 = (vf. ) G by 
setting 

4, P (t) = u v Kp {t) -e if t G i? > 
= «*,,(*) + £ if 

and setting v £ hr (t) = u" ir (t) in all other cases. Notice that, if e = min{u£ pi , 4« — •} then 
these new departure rates are still admissible. By (4.37) and (4.38), this construction yields 



h:r 



E / CW-krW-^))* = e I ^Z\t)dt-e I *%>(t)dt < -2e6, (4.39) 



providing a contradiction with (4.32). 



CASE 2: u? ■ = 4k. If this equality holds, consider the index 

j* = max {i < j ; u£ ^ < 4k} . 



Notice that tj* > — T. Indeed, by construction t 2 > -T . If t v r < -T, by (4.28) this would 
imply 

u v k q (t) = 4k for all t G [f£. ,t£] 5 [-T, -T ] , 

ju u Kq (t)dt > 4/e(tjf-i£) > 4k(T-T )> G, 
reaching a contradiction. We consider two subcases. 

CASE 2a: = j — 1. In this case, since it is not restrictive to assume — < |, we have 
Jj'_ 1 = [ij_2i ^f-i] c — 8, t 2 + 5]. We can thus derive a contradiction as in CASE 1, simply 
replacing j by j — 1. 



CASE 2b: j* < j - 2. Observe that, for all s x < s 2 , 

1 fsa 

7fc.a(*2) - r fc „(si) > 



± max J si 



In particular, for any s\ G ij* and S2 G /J we have 



Tk,qM -l~k, q (si) > 

This yields the estimate 



* 2 



- 1 max Jsi 



4k [tj-i > 



4k(s 2 - gi) 



""mm 



(4.40) 



(4.41) 



(4.42) 
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On the other hand, we have 

<Pk(s 2 ) ~ ipk(si) > - y'maxisi ~ s x ) . 
Recalling the definition of the constant k in (4.27), from ( |4.42 )-(4.43) we obtain 

'4/c0' „•„ A , I 



3R 



mm / 

rmax 



(S2 - Si) 



(S2 - Si) 



(4.43) 



(4.44) 



for all si € ij» and s 2 G ij . 

We now choose the departure rates v £ = (f| p ) G £/ w by setting 

= " e if * e IJ , 

vl, q (t) = if t€ij., 

and setting v e hr {t) = u^ r (t) in all other cases. Notice that, if e 
then these new departure rates are still admissible. 



(4.45) 



Using (4.44) with si = t, S2 = * + P- — t v - t we compute 

cT 



E 



T 



h,r 



(t) ■ lvl r {t) - u v h Jt)\ dt 



1% 

3* 



(a"), 

k,q 



(ty-ts.) < o. 



Once again we reached a contradiction with (4.32), completing the proof. 



(4.46) 



□ 



Remark 4.4. The above proof is based on a topological method, and does not yield any 
information about the uniqueness of the Nash equilibrium. Another important question is the 
dynamic stability of this equilibrium solution. This issue was investigated numerically in [3]. 
In the case of a single group of drivers traveling on a single road, the uniqueness of the Nash 
equilibrium solution was proved in [2], but the stability issue remains unresolved. 
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